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Abstract As an observational case study, we consider the origin of a prominent 
poleward surge of leading polarity, visible in the magnetic butterfly diagram 
during Solar Cycle 24. A new technique is developed for assimilating individual 
regions of strong magnetic flux into a surface flux transport model. By isolating 
the contribution of each of these regions, the model shows the surge to originate 
primarily in a single high-latitude activity group consisting of a bipolar active re¬ 
gion present in Carrington Rotations 2104-05 (November 2010-January 2011) 
and a multipolar active region in Rotations 2107-08 (February - April 2011). 
This group had a strong axial dipole moment opposed to Joy’s law. On the 
other hand, the modelling suggests that the transient influence of this group on 
the butterfly diagram will not be matched by a large long-term contribution to 
the polar held, because of its location at high latitude. This is in accordance 
with previous hux transport models. 

Keywords: Active Regions, Magnetic Fields; Magnetic helds. Photosphere; 
Solar Cycle, Models 


1. Introduction 


The gradual transport of magnetic hux from decaying active regions on the 
Sun’s surface toward the polar re gions is essential for explaining the observed 
pattern of surface magnetic held ( Babcockl. Il96ll: Ibeightonl. Il964 Wang, Nash , 
and Sheelev. Il989bl: Mackav and Yeated . l2012l : IPetrie. Petrovav. and Schattenl . 
I20l4ljiang et al . 2014ll . However, the net magnetic hux transported to the poles 
plays an even more signihcant role in the popular Babcock-Leighton solar cycle 
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model ( Babcock . 1961 : Leighton . 1969l l. In this picture, the average latitudinal 
offse t between positive and negative active region polarities (known as Joy’s 


law: Hale et ~al\ . 1919l l leads to an imbalance between the overall poleward 


transport of positive and negative polarities in each hemisphere. This forms a net 
poloidal magnetic field, and differential rotation of this field in the underlying 
convection zone is believed to produce the azimuthal magnetic field driving the 
subsequent solar cycle of sunspot emergence. This viewpoint is supported by 
observed correlations between the polar field at solar minimum and the strength 
of the subsequent solar cycle ('e.o. iMunoz-Jaramillo et aH 2013l l. Unfortunately, 
our ability to predict future solar cycles remains limited by our ability to predict 
differences in the poloidal field production from one cycle to the next. This 
motivates continuing study of surface-flux transport on the Sun . _ 

It has been recognised since the work of Howard and Labontel (1981) that the 
transport of magnetic flux from decaying active regions to the poles takes place 
not smoothly, but through episodic “surges”. This is seen by inspecting the mid¬ 
latitudes in Figure [T^, which shows the so-called butterfly diagram of longitude- 
averaged radial magnetic field on the photosphere, for Solar Cycle 24. The data 
used here are radial-component [Bj.] synoptic maps from the Global Oscillation 
Network Group (GONG, http://gong.nso.edu/data/magmap/). In particular, we 
focus in this article on a particular poleward surge of negative polarity in the 
northern hemisphere, labelled A in Figure[I}i. This surge had the leading polarity 
for this hemisphere, meaning the polarity of (the majority of) equatorward spots 
in bipolar groups. 

In Section [2l we describe the data-driven surface-flux transport model used 
to analyse the origin of surge A. The longitude-averaged Br for the simulation 
is shown in Figure [Ud. In fact, polar surges wer e important in the historical 
development of the surface-flux transport model ( Sheelevl . l2005l l. because they 
represent strong evidence for the presence o f a poleward meridional flow. This 
was suggested bviHoward and Labontel (|l98in . and demonstrated conclusively by 


Wang. Nash, and Sheelevl l|1989al ). who used surface-flux transport simulations 


to model Solar Cycle 21. Their model reproduced surges only when they included 
a poleward meridional flow in addition to supergranular diffusion. 

The timing of the poleward surges is dependent both on the flux-transport 
parameters and, of course, on fluctuations in the underlying pattern of flux 
emergence. For example, using both Ki tt Pea k and Mount Wilson Observa¬ 
tory magnetograms from 1974 to 2012, Petriel ( 2012h showed that the yearly 
averaged magnetic-field strength around 50° latitude closely tracks a measure 
of total active region flux, provided that the latter is w eighted by the latitude- 
differe nce between centroids of pos itive and negative flux. IZolotova and PonvavinI 
( 2 OI 2 I I developed an earlier idea of Gnevvshev ( 1938h that the 11-year solar cycle 
consists of two or more superimposed activity impulses, different in each hemi¬ 
sphere. They suggested that these impulses produce poleward surges of trailing 
polarity in each hemisphere, while allowing waves of leading polarity (such as 
our surge A) to occur in the gaps between impulses. The impulses themselves 
ma y correspond to the complexes of activity studied by Gaizanskas et al\ 1 198,81 1 
and Gaizausl^ ( 2008[l . which comprise multiple emerging active regions at the 
same location over three to six solar rotations. Observations suggest that these 
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Figure 1. Longitude-averaged Br as a function of time and latitude, in (a) the sequence 
of observed GONG synoptic maps and (b) the flux-transport simulation (Section [2|| for the 
period GR 2074-GR 2154 (August 2008-August 2014). Both figures are saturated at ±3 G. 
The leading-polarity (negative) surge A is indicated in panel (a). 


impulses are also evident in coronal extreme-ultraviolet and soft X-ray emission 
( Benevolenskaval . 2003h . In Section [31 we show that the leading-polarity surge 
A originated primarily from a single active region during Carrington Rotation 
2107. It is better interpreted as falling between two successive impulses, rather 
than itself being associated with a long-lived complex of activity. 

A common property of poleward surges is that t hey tend to occur in p airs of 
opposite polarity. This is clearly seen in Figure 8 of lUlrich and Tranl ( 2013 ). who 
plot dBr/dt from Mount Wilson Observatory measurements and see a pattern 
that they call “ripples”. The successive arrival of alternating polarity surges 
causes significant short-term oscillations in the polar field. When the contribu¬ 
tion of an individual active region is isolated (as we will illustrate in Section|3]), it 
tends to contribute to a poleward surge of each polarity, typically offset in time. 
Later in the cygls^T^ggg^g^iibuligaiJaigglijgaSggLgnSjASglhs^Wang, Nash, 
and Sheelev Jl989ah . As we will discuss further in SectionHl the long-term impact 
on the polar field then depends on the imbalance between these contributions, 
which is ultimately equivalent to the imbalance of net magnetic flux from the 
region in each hemisphere, which in turn results from cross-hemispheric flux 
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transp ort and magnetic cancellation iGiovanellil. 19851: Durrant. Turner, and 
Wilson, 2004 ICameron et all . 201 , Sl l. 


2. Data-driven Flux Transport Model 


We model the global photospheric radial mag netic field, \Br(0, (j), t)], during Solar 
Cycle 24 using a surface- flux transport model ( Sheelev , 2005 ; Mackav and Yeat^ . 


2012 : Jiang et'^ . l20l4 ). In particular, we write 


in spherical polar coordinates (r, 9, ^), and evolve the vector potential [{Ag, A^)] 
numerically on r = Rq with a finite-difference method through the equations 


dAe 

dt 

dA^ 

dt 


T) f)fi 

a;(6»)i?o sinOBr - ——+ -S'e(6», </>,t), 
/tQSine/ 0(j) 

T) FiFi 

— UgBr + 0 - 


( 2 ) 

(3) 


The equations are solved in the Carrington frame, in which we take the angular 
velocity of differential rotation to be 

u}{9) = 0.521 — 2.396 cos^ 9 — 1.787cos^ 9 degday”^. (4) 

The meridional circulation uses the Schiissler and Baumann! (2006) form 

= “0 sin(2A) exp (tt — 2|A|), (5) 


where A = n/2 — 9 is latitude. In this article, we choose uq constant in time 
to give a peak flow speed of 11 ms“^. We also choose a constant supergranular 
diffusivity of 77 = 500 km^ s“^. These values were chosen to optimize agree ment 
between the simulated and observed butterfly diagrams ( cf. lYeatesl . l2014l l . For 
simplicity, we do not include any additional flux-decay term. The terms Sg{9, </>, t) 
and S(i,{9, (j), t) in Equations dH) and ([2]) represent the emergence of new magnetic 
flux from beneath the photosphere. These terms are treated in a novel way, as 
described in Section HD below. 

For this article, we simulate the period between Carrington Rotations 2074 
and 2154 (namely 30th August 2008 and 17th September 2014). The observa¬ 
tional input data used to generate both the initial Bj. map and the source terms 
80 ( 9 , (pA), Sci,{9,(l),t) are radial-component synoptic maps from GONG. These 
were chosen because they provide continuous coverage with consistent calibration 
over our period of interest, at a suitable resolution (360 pixels in (p and 180 pixels 
in sin A). For simplicity, we use a single synoptic map per Carrington rotation, 
with the most recent update at the end of the rotation. This time-resolution 
of input data is sufficient for present purposes. As already mentioned above. 
Figure [T] shows the resulting “butterfly diagram” of longitude-averaged Br on 
the photosphere, for both the observed magnetograms and the simulation. 
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2.1. Assimilation Technique 


Our treatment of the emerging flux terms Ss{0, (j), t) and S^{9, (j), t) in Equations 
© and ([3]) differs from previous models. Typic ally, new flux is added either in the 
form of discrete bipolar i nagnetic regions (e.q. Sheelev. DeVore. and Bori^ . 1985 : 
Mackav and Lockwoodl. 12002 JSchriiver 4 _DeJRpsaj,^nd_Titl^ 20021: Biiu mann, 


Sch mitt, and Schiissler.l200^ Yeates. Mackav. and van Ballegooiienl.l2007l: Jiang 

et al. j2010all or bv assimilating observed magnetograms on the near-s ide of the 
Sun ( Worden and Harvev . 20^: ISchriiver. De Rosa, and Titlel. l2003l: Durrant, 


Turner, and Wilson, 2004; Arge et ai . 201fll : Upton and Hathaway! . ^Idh . The 
latter leads to more realistic reproduction of the observed photospheric held, 
but makes it difficult to isolate the contribution of individual emerging regions, 
as is our goal in this article. On the other hand, previous models using bipolar 
magnetic regions have assumed an idealized form for the newly emerging regions 
that is limited in accuracy for complex active regions. As a result, we employ 
here a new technique for assimilating individual strong-hux regions into the 
simulation. The aim is to separate individual new regions, while maintaining 
their detailed structure from the observed magnetogram. 

The assimilation algorithm is described in detail in Appendix [A] To illustrate. 
Figure [2] shows the strong-hux regions identihed by our numerical code for CR 
2132. As can be seen, the algorithm identihes both isolated bipolar magnetic 
regions and multipolar complexes of activity. The regions are incorporated one- 
by-one into the simulation on their day of central-meridian crossing. 


■■= 0 


-50 



Figure 2. Strong-flux regions identifled by the automated technique for the GONG synoptic 
map of GR 2132. Region outlines are shown in different colours, superimposed on the observed 
synoptic map of Br{6,4)) (saturated at ±30G). 

There are two major advantages of this procedure over previous methods 
for assimilation of purely bipolar regions. Firstly, it is fully automated, giving 
consistency and objectivity in selection of new hux regions, as well as saving a 
considerable amount of time. Secondly, the new strong-hux regions replace the 
pre-existing simulated hux within those regions, rather than being superimposed. 
This is a major advantage because it allows the model to incorporate subsequent 
hux emergence observed over multiple solar rotatio ns within exis ting activity 
complexes. In the bipole-based emergence model {e.g. Yeatesl . 2014l l. this is very 
difficult to model realistically. Yet the structure of these activity groups plays a 
major role in the global distribution of magnetic hux on the Sun. 
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2.2. Grouping of New Regions 

For the analysis in this article, we make the additional step of grouping together 
new regions when they emerge at the same location in space, separated in time 
by at most two Carrington rotations. This is necessary because the assimilation 
algorithm often modifies an already emerged region, if it still has strong flux in 
a subsequent rotation. Each group should therefore be considered together when 
assessing its contribution to a polar surge or the polar field. Of the 204 strong- 
flux regions identified between CR 2074 and CR 2116, we identify 96 separate 
groups. The groups range in size from 1 to 15 regions, with a mean size of 2.125 
regions, but large groups are rare and the median size is only 1 region. 


3. Origin of Poleward Surge A 

To determine which of the emerged-flux groups produced polar surge A, we 
repeated the flux-transport simulation omitting each of the groups in turn. We 
then quantify the effect of each group by cross-correlating the butterfly diagram 
from each run against that from the complete run, in a window from CR 2104 
to CR 2117 and A = 35° to A = 70°. This window is chosen to pick out surge A. 

For example. Figure [3] shows the butterfly diagram when group 51 is omitted, 
compared to that from the complete run. The window cross-correlation between 
the complete run and the observed butterfly diagram (Figured^) is ro = 0.864. 
This residual error must result from processes that are not fully accounted for in 
our simple flux-transport model. By contrast, the run with group 51 omitted has 
a lower cross-correlation, r = 0.682, indicating that group 51 has a significant 
effect on surge A. 

In fact, the correlation analysis shows that group 51 is primarily responsible 
for surge A. To see this. Figure H] shows the cross-correlation coefficient [r] for 
all groups (up to and including CR 2116). The colour and size of each symbol 
denote the average axial dipole moment of each group, to be discussed below. 
Focusing for now on the r-values, it is seen that omitting the majority of groups 
has no effect on the correlation in this window. Only groups 22, 40, 47, 51, and 
52 cause r to differ from tq by more than 0.02, and only group 51 causes a 
difference of more than 0.05. For reference, the NOAA active-region numbers 
associated with these five groups are given in Table [TJ 

Before considering group 51 in more detail, let us remark on the other five 
groups affecting our correlation window. Groups 22, 40, and 47 each consist of 
only a single bipolar region, tilted in accordance with Joy’s law. These regions 
contribute to the cross-correlation primarily by strengthening the positive flux 
region in the top left of the correlation window, although their leading polarity 
flux does contribute to the negative surge A. Group 52 is notable in that omitting 
it actually increases r [i.e. improves the match with the observed butterfly 
diagram). This group is a large cluster of multiple bipolar regions, persisting 
over seven rotations (CR 2104 to CR 2111). It is responsible for much of the 
positive surge in the bottom-right of the correlation window. Compared to the 
observations, this surge is too strong; thus omitting group 52 actually improves 
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Figure 3. Effect of omitting group 51 from the flux-transport simulation, shown in longi¬ 
tude-averaged Br- The top panel shows the run with group 51 omitted, while the bottom 
panel shows the complete run. The cross-correlations r and ro are computed in the indicated 
window, against the observed GONG butterfly diagram. 


Table 1. NOAA active-region numbers for the groups affecting the correlation 
window. 


Group 

Carrington rotations 

NOAA regions 

22 

GR 2093 

11048 

40 

CR 2100 

11097 

47 

CR 2103 

11120 

51^ 

CR 2104-5, CR 2107-8 

11131, 11140, 11164, 11180 

52 

CR 2104, CR 2106-11 

11135, 11157, 11169, 11185-6, 11189, 
11203-5, 11228-30, 11242 


^Illustrated in Figure 5. 


the correlation. This suggests that the processes in our flux-transport model may 
be overestimating the poleward spreading from large complexes such as group 
52, per haps because we neglect the effect of systematic inflows towards activ e 
regions ( De Rosa and Schriiverl . [200^ : lBasu and Antial . 201fll : Jiang et al. . 2010tjl . 

The majority of the negative flux in the poleward surge clearly arises from 
group 51. This group consists of four individual strong-flux regions, determined 
from the synoptic maps for CR 2104, CR 2105, CR 2107, and CR 2108. These are 
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Figure 4. Cross-correlation between the observed and simulated butterfly diagrams. Each 
triangle corresponds to a simulation with a different group of regions omitted, located at 
the time of insertion of the first region in that group. The colour and size of the triangle are 
proportional to the (average) axial dipole moment [G] among that group of inserted regions, as 
indicated in the colourbar, with upward triangles for positive values and downward for negative. 
Dashed-vertical lines show the extent of the correlation window, while the dashed-horizontal 
line shows the cross-correlation [ro] when all groups are included. 

highlighted in the synoptic maps in Figure [SJ above the simulated field (at the 
end of the corresponding Carrington rotation). It is seen that group 51 consists 
of i) a bipolar active region appearing first in CR 2104, and still having strong 
flux in CR 2105, and ii) a multipolar region appearing first in CR 2107 and still 
visible in CR 2108. Both of these constituent regions appear to have a large 
tilt angle, which is understood as the angle between their positive and negative 
centroids with respect to the East-West line, although they are oppositely tilted. 
Since the later region is multipolar in nature, a better way to quantify this is 
the axial dipole moment of the region, defined as 



( 6 ) 


The four strong-flux regions comprising group 51 have 6i,o of 0.0120 G, 0.0182 G, 
—0.0546 G, and —0.0205 G respectively. Hence the group average is negative 
(—0.0112 G), as shown in Figure 01 We argue that this relatively large and 
predominantly negative dipole moment of the later region in group 51 (opposite 
to the majority positive sign from Joy’s law) is primarily responsible for surge A 
in the butterfly diagram. The effect is strengthened by the relatively high latitude 
at which this region emerges. However, there is also a smaller contribution from 
the leading polarity of the earlier bipolar region in the group, and to a lesser 
extent from groups 22, 40, and 47. 
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Figure 5. The four strong-flux regions comprising group 51, shown outlined in the GONG 
synoptic magnetograms (top row) and simulated Br (bottom row). The simulated Br is shown 
at the end of the corresponding Garrington Rotation. All images are saturated at ±15 G, yellow 
positive and blue negative. 


4. Contribution to the Polar Field 


It would be erroneous to suppose that, merely because group 51 causes a partic¬ 
ularly prominent opposite-polarity polar surge in the butterfly diagram, it must 
make a particularly strong contribution to the net polar field. In fact, since group 
51 is located at rather high latitude, the flux transport model predicts that its net 
polar field contribution at the end of Cycle 24 will be weak. The strong poleward 
surge A seen in the butterfly diagram is merely a transient phenomenon. This is 
well known from previous studies of surfac e-flux transport (IWang and Sheelevl . 
1991 : Mackav. Priest, and Lockwoodl . [2002h . but it is interesting to illustrate this 


for our specific region 51. 

To see the relative effect of the different emerging groups on the polar field, 
we show in Figure [S] the change in the north polar flux 


^Npit) = Rl[ dn (7) 

Je<20° 

as each group is omitted from the simulation. It is seen that the largest changes 
of polar flux come from a relatively small number of emerging groups, and many 
of these are the groups with largest |5i,o|- (The others are groups with many 
constituent regions, but weak average dipole moment.) Additionally, the sign of 
the difference in Figure[6]is seen to correlate with the sign of bi^Q for each group. 
The majority of groups have 6i o > 0 (in accordance with Joy’s law for Cycle 
24), and therefore contribute positively to d’Np. Group 51 is conspicuous in its 
negative and one of few groups to make a net negative contribution to $np. 

A crucial observation in Figure [5] is that the polar-field contribution from 
many groups - including group 51 - appears to initially peak, then to decay in 
time. This is entirely typical: the initial peak is caused by the negative surge 
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Figure 6. Change in north-polar flux [^np(^)] compared to the full simulation, in the runs 
where each emerging group is omitted. The curves are coloured according to the average axial 
dipole moment [6i,o] of the group [G]. The curve for group 51 is labelled. 


reaching the pole first, and then the decay is caused by the positive surge from 
the same region reaching the pole at a later time. 

It is cleanest to illustrate this effect by simulating the evolution of a single 
regioninjsolationj^^ithoirttiie^urroundin^^eld^X£i 2 ^^^®'Ckay, Priest, and Lock- 
wood, 2002h . We choose the strong-flux region emerging in CR 2107 (i.e. in the 
third column of Figure [51). Figure [7] shows three butterfly diagrams obtained in 
such single region simulations, when the complex is placed i) at its observed 
latitude A = Aq = 24°N, ii) at A = Ao/2, and iii) at A = 0. For the same three 
runs. Figure [5] shows both the polar flux [$np( 0] the axial dipole [l>i,o(^)] 
as functions oFtime^^heJattei_can,becom 2 ared^o^dgure 6 ,rf Jiang, Cameron, 
and Schiissler (l2014l i and Figure 4 of Wang and SheelevI (|l99ll l. 

From Figures [7] and we see the following behaviour: 

i) The polar flux does not (in general) decay to zero, but reaches an asymptotic 
steady state after about six years, which is symmetric between the two 
hemispheres. (This results from a bala nce between equatorward diffusiv e 
transport and poleward meridional flow: lDeVore. Boris, and Sheelevl ■ ll984^ . 


ii) Placing the emerging region nearer the Equator leads to a larger asymptotic 
value of the polar field. (The exact latitudinal dependence depends on the 
chosen meridional-flow profile and supergranular diffusivity). 

iii) For the latitudes shown, the axial dipole moment does not show the turnover 
seen in the polar field for Aq and Ao/2. This is because it does not have 
to wait for flux to be transported to the pole to register the presence of 
the region. For the same reason, it tends to converge more rapidly to its 
asymptotic value than does the polar flux, although the asymptotic values 
of the two quantities are proportional. 

It is clear that, despite its significant contribution to the transient poleward 
surge A, group 51 is likely to make a relatively small contribution to the polar 


SOLA: ms-v2.tex; 18 February 2015; 1:20; p. 10 





















Source of a Prominent Poleward Surge 


D) 

CD 

■D 


D) 

a> 

■D 



Years 


Figure 7. Longitude-averaged Br in the three simulations of an identical region emerged at 
different latitudes. (Yellow shows positive, blue negative). 


field at the end of Cycle 24 owing to its high latitude of emergence. The key 
point is that a decaying active region will produce a net (asymptotic) polar 
field only i f it leads to a net imbalance in magnetic flux in e ach hemisphere 
( Giovaneiiil . 19851 : IPeVorel . Il98^ : Wang. Lean, and Sheelev . [2002h . This can only 
happen if net flux from the region either emerges on both sides of the Equator, 
or is transported across it. It follows that those active regions located near the 
Equator - and in particular straddling the Equator - are the mo st important for 
prod ucing fluctuations in the polar-field production process (c/. ICameron et al 
1201,811 . 


5. Conclusions 


In summary, the activity group responsible for poleward surge A is a striking 
example of the statistically fluctuating nature of active-region flux emergence on 
the Sun. This group has a large axial dipole moment, opposite in sign to Joy’s 
law. It is the primary source of the poleward surge A that is clearly seen in the 
magnetic butterfly diagram. 

Erom the viewpoint of the Babcock-Leighton model for the solar cycle, the 
importance of an individual emerging region lies in its net contribution to the 
polar field. Our modelling shows that this particular activity group, whilst caus¬ 
ing a significant transient perturbation to the polar field, will not lead to a large 
long-term contribution. This is because it is located far from the Equator, and 
it is the cross-equatorial fl ux transport that controls the asymptotic polar field 
(as was first recognised bv iGiovanelli . 1985). On the other hand, the transient 
perturbation does last for about six years. This suggests that higher-latitude 
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Figure 8 . North polar flux [<t’jvp(4)] and axial dipole moment [ 61,0 (4)] as a function of time 
for the three simulations of an identical region emerged at different latitudes. 


regions emerging later in the cycle may still be significant contributors to polar 
field measurements taken at a fixed time late in the cycle. It is also not clear that 
the best surface indicator of underl ying poloidal magnetic field in the convection 
zone will be the polar field alone (jWang and Sheelevl . Il99lh . The sub-surface 
decay of active-region flux at lower latitudes may well be an important f actor in 
repairing the azimuthaI_^eM_beltsj|vaji_^aJle 20 oijen^_and_^laefe^200^Yeates 
and Muhoz-Jaramillo, 2013ll . In this case, the role of high-latitude active regions 
in the Babcock-Leighton process may be more significant than their polar-field 
contribution might suggest. 

Finally, we have developed a new technique for data-assimilation in the flux 
transport model. This demonstrates the feasibility of a middle ground between 
inserting idealised bipolar regions and assimilating full magnetograms. In the 
future this could be combined with more sophisticated automated techniques 
for detecting_j,ndJnserting,^ctive,j;e£ions,^d,Jnghei_^ime_cadence_[ei 2 ^IIiggins 


et al.. 
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Appendix 

A. Assimilation Algorithm 

To assimilate new magnetogram flux, we first process each individual synoptic 
magnetogram 5 ( 0 , (j)) to determine a list of strong-flux regions, using the following 
procedure: 

i) Correct any overall net flux imbalance in g by subtracting an equal amount 
from each pixel. 

ii) Take the absolute value to obtain a new map = \g{9,(f))\. 

hi) Smooth / with a Gaussian filter (standard deviation cr) to obtain /. The 
purpose of this step is primarily to merge together the positive and negative 
polarities of individual bipolar active regions. 

iv) Identify connected regions [i?i] where f > Bq- These are the new regions 
that will be inserted into the flux-transport model. (Note that these regions 
need not be simply connected and could have holes.) 

Next, each strong flux region [i?i] is incorporated into the flux-transport model 
as follows: 

i) Determine the centroid of Ri in Carrington longitude, and choose to insert 
it on the corresponding day when this longitude crosses central meridian. 

ii) Determine the net flux imbalance B^dfl, where Br{0,(j)) is the 

pre-existing simulated radial field. 

hi) Modify the observed field so that the net flux in Ri matches by 

subtracting an equal amount from each pixel. This prevents us introducing 
a flux imbalance in the insertion step. 

iv) Replace the simulated Br in Ri with the observed held 

v) Recompute the global vector potential [Ag, A^] to match the modihed Br- 

The process is designed to be fully automated, once the controlling parameters 
(T and Bq have been determined. For the GONG synoptic maps, we have found 
a = 3 pixels and Bq = 15 G to work well. The smoothing width a should be 
large enough to merge the two opposite polarities of activity complexes, while 
not merging too many neighbouring complexes. The threshold Bq needs to be 
large enough to include all signihcant new flux, but be above the “noise” level of 
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ephemeral regions (in the smoothed magnetogram), so that only flux at active 
latitudes is assimilated. 

One must be aware that, in order to preserve the field distribution, the net 
flux in each identified region Ri should be close to that existing beforehand in the 
simulation. For example, a region emerging in an empty area of the solar surface 
should have zero net flux. Otherwise, the flux correction will noticeably distort 
the observed field. This practical consideration requires cr to be sufficiently large 
and Bq to be sufficiently small. 
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